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ABSTRACT 

We give an overview of our recent integral- field-unit spectroscopy of luminous extended emission-line re- 
gions (EELRs) around low-redshift quasars, including new observations of 5 fields. Previous work has shown 
that the most luminous EELRs are found almost exclusively around steep- spectrum radio-loud quasars, with 
apparently disordered global velocity fields, and little, if any, morphological correlation with either the host- 
galaxy or the radio structure. Our new observations confirm and expand these results. The EELRs often show 
some clouds with velocities exceeding 500 km s~\ ranging up to 1100 km s~\ but the velocity dispersions, 
with few exceptions, are in the 30-100 km s~^ range. Emission-line ratios show that the EELRs are clearly 
photoionized by the quasars. Masses of the EELRs range up to > 10^^ Mq. Essentially all of the EELRs show 
relatively low metallicities, and they are associated with quasars that, in contrast to most, show similarly low 
metallicities in their broad-line regions. The two objects in our sample that do not have classical double-lobed 
radio morphologies (3C48, with a compact- steep- spectrum source; Mrkl014, radio-quiet, but with a weak 
compact- steep- spectrum source) are the only ones that appear to have recent star formation. While some of 
the less-luminous EELRs may have other origins, the most likely explanation for the ones in our sample is 
that they are examples of gas swept out of the host galaxy by a large-solid-angle blast wave accompanying the 
production of the radio jets. The triggering of the quasar activity is almost certainly the result of the merger of 
a gas-rich galaxy with a massive, gas-poor galaxy hosting the supermassive black hole. 
Subject headings: quasars: individuals(3C48, Mrk 1014, Ton 616, Ton 202, PKS 2251+11, 3C 249.1, 4C 37.43, 

3C 79) — quasars: emission lines — galaxies: evolution — galaxies: ISM — galaxies: 

abundances 



" ; 1. INTRODUCTION 

\ As soon as quasars were discovered, it was noticed that 
some of them showed detectable nebulosity surrounding 
their star-like nuclei (e.g., Matthews & Sandage 1963). In 
an attempt to understand the nature of the fuzz around 
3C48 (Matthews & Sandage 1963; Sandage & Miller 1966), 
Wampler et al. (1975) serendipitously discovered extended 
QQ \ emission-line gas in a spectrum taken at a region r^A" north 
of the quasar, although this "fuzz" they intended to look into 
1^ turned out to be mostly the unusually large and disturbed 
. 1^ host galaxy of 3C48 (Boroson & Oke 1982). Shortly after- 
s' , ward similarly extended nebulae were identified in off-nuclear 
^ ' spectra of two other quasars: 4C 37.43 (Stockton 1976) and 
_C3 ■ 3C 249. 1 (Richstone & Oke 1977). In spite of the limited spa- 
tial sampling of these observations, it was already clear that 
such nebulae extend well beyond the classical narrow-line re- 
gion (NLR) which is typically confined within ~ 1 kpc from 
the nucleus. These "extended emission-line regions" (EELRs; 
following Stockton & MacKenty 1983), are clearly physically 
associated with the quasars: the ionized gas shows almost the 
same redshift as the nearby quasar, and the spectrum of the gas 
suggests the quasar as a more likely photoionization source 
than massive young stars. 

Slit spectroscopy does not allow an accurate determination 

^ Based in part on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the 
Particle Physics and Astronomy Research Council (United Kingdom), the 
National Research Council (Canada), CONIC YT (Chile), the Austrahan Re- 
search Council (AustraHa), CNPq (Brazil) and CONICET (Argentina). Gem- 
ini Program ID: GN-2007A-Q-43 and GN-2007B-Q-12. 

^ Current address: Department of Astronomy, California Institute of Tech- 
nology, MS 105-24, Pasadena, CA 91125; fu@astro.caltech.edu 



of the morphology of the extended nebulae. As the redshift 
of 3C 249.1 conveniently places the [Oiii] A5007 line into 
the commonly available Ha filter, its EELR was the first to be 
imaged through a narrow-band filter centered on this strongest 
optical nebulae line (Stockton & MacKenty 1983). The gas 
was found to form two extensive "tidal-tails" to the NW and 
SE of the nucleus. As more narrow-band filters were ac- 
quired, Stockton & MacKenty continued the imaging survey 
on a sample of 58 quasars at z < 0.5. By the end of the survey, 
47 quasars were imaged in [O III], and extended emission was 
detected in 15 of them (Stockton & MacKenty 1987). These 
images reveal that the morphologies of EELRs are often fas- 
cinatingly complex and clearly betray a non-equilibrium situ- 
ation showing knots and filaments that sometimes extend for 
tens of kpc. Structures reminiscent of tidal tails such as the 
EELR of 3C 249.1 turn out to be rare — Ton 202 is the only 
other similar case known. Surprisingly, the distribution of the 
ionized gas generally bears no close morphological relation- 
ships either with the host galaxies or with the extended radio 
structures, if present. 

The origin and nature of these impressive and often beau- 
tiful structures have proved to be elusive. In this paper, we 
present new imaging spectroscopic observations of 5 quasars 
with luminous EELRs, discussing these in the context of our 
more detailed investigations of 3 additional EELRs (Fu & 
Stockton 2006, 2007b, 2008) and the surprising correlation 
found between quasar broad-line metallicity and the presence 
of luminous EELRs (Fu & Stockton 2007a). Section 2 gives 
an outline of previous work on EELRs, § 3 describes our ob- 
servations and data-reduction procedures, § 4 gives our results 
for the 5 quasars with new observations, § 5 discusses the ion- 
ization mechanism and the metallicities for these EELRs, and 
§ 6 gives an over-all summary, including a general scenario 
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that we believe can account for the observed properties of 
luminous EELRs. Throughout the paper, we assume a con- 
cordance cosmological model with = 70 km s~^ Mpc~\ 
ftm = 0.3, and = 0.7. Unless stated otherwise, solar abun- 
dances are defined by the results in Anders & Grevesse (1989) 
[12+log(O/H)0 = 8.93 and Z© = 0.02]. Even though these 
are not the most recent determinations of solar abundances, 
they have the advantage that they allow easy comparison with 
many previous studies that have used this metallicity scale. 

2. BACKGROUND 

2.1. Correlations between Extended Emission and Quasar 
Properties 

Boroson & Oke (1984) and Boroson et al. (1985) carried 
out an important spectroscopic survey to obtain off-nuclear 
spectra of 9 radio-loud quasars and 3 QSOs^. In addition to 
finding more quasars with extended emission, an important 
discovery of this program was the strong correlation between 
luminous optical extended emission and radio morphology 
(or radio spectral index). Boroson et al. found that all of 
the six quasars with strong extended emission lines are steep- 
spectrum radio sources, and the other six objects'^ where such 
extended emission is weak or absent are either radio-quiet or 
flat- spectrum radio-loud sources. There is a strong correlation 
between radio spectral index and radio morphology, since ex- 
tended radio structures such as radio lobes are dominated by 
optically thin synchrotron radiation, which produces a steep 
power-law spectrum, while unresolved radio cores gener- 
ally shows a flat spectrum due to synchrotron self- absorption 
(or thermal bremsstrahlung radiation from a disk wind, see 
Blundell & Kuncic 2007). This extended emission-spectral 
index correlation can thus be translated into an extended 
emission-radio morphology correlation — luminous extended 
emission was found exclusively in quasars having extended 
radio structures, most of which show Fanaroff-Riley type II 
(FR II) twin-jet morphologies (Fanaroff & Riley 1974). 

This extended emission-radio morphology correlation was 
confirmed by the [Olll] narrow-band imaging survey of 
Stockton & MacKenty (1987), but it evolved into a more 
complicated picture as a result of their 4 x larger sample and 
their seamless spatial sampling. The authors found that (1) 
not all steep- spectrum {a^ > 0.5) lobe-dominated quasars 
are surrounded by luminous EELRs — only 10 out of the 
26 (38 ± 12%) have EELRs with Le[oiii] > ^cut = 5 x 10"^^ 
erg s~^ (hereafter EELR quasars), where Le[oiii] is the total 
[Olll] A5007 luminosity within an annulus of inner radius 
11.2 kpc and outer radius 44.9 kpc centered on the nucleus 
and Lcut is determined by the highest Le[oiii] upper limits 
of the objects where extended emission is not detected (see 
Fig. la)', and (2) none of the seven flat-spectrum (a^ < 0.5) 
core-dominated quasars and only one of the 14 (7 ± 7%) op- 
tically selected QSOs show an EELR with Le[oiii] > ^cut- 
To illustrate this correlation more clearly, we have created a 
new figure (Fig. la) based on the original data of Stockton & 
MacKenty (1987). 

The sample observed by Stockton & MacKenty (1987) in- 
cluded almost 2/3 of all optically luminous QSOs known at 

^ Here we consider only their high-luminosity sample to be consistent with 
the sample of Stockton & MacKenty (1987) which we will discuss in the 
following paragraphs. 

^ Contrary to these original papers, 3C 48, a compact steep-spectrum radio 
source with a powerful radio jet extending 0^/1 to the north, is in fact sur- 
rounded by a luminous EELR, as was later shown by Stockton & MacKenty 
(1987). 



the time with z < 0.45 and 6 > -25°. It is still nearly complete 
for radio-loud quasars within these limits but clearly incom- 
plete for radio-quiet QSOs, which typically comprise about 
90% of the total population. Nevertheless, the radio-quiet ob- 
jects in the sample are likely typical of the population of lu- 
minous, low-redshift, UV-selected QSOs. 

Boroson et al. (1985) also noticed that the presence of 
strong extended emission was frequently accompanied by 
broader and bumpier nuclear B aimer lines, weaker nuclear 
Fe II emission, and stronger nuclear narrow lines. The first 
two correlations are often regarded as by-products of the ex- 
tended emission-radio morphology correlation, as there exist 
well-established correlations between radio morphology and 
nuclear broad-line profile and Fe II emission (Miley & Miller 
1979; Steiner 1981). The last correlation seems to be an in- 
trinsic one, and it was confirmed by Stockton & MacKenty 
(1987), who found that strong extended [Olll] emission is 
present only for objects with strong nuclear narrow-line emis- 
sion. Specifically, all of the 11 EELR quasars in Stockton 
& MacKenty (1987) have nuclear [Oiii] A5007 luminosity 
(^N[oiii]) of at least 6.5 x 10"^^ erg s~\ which is clearly seen 
in Fig. lb. 

The nuclear [O III] A5007 luminosity has frequently been 
used as a proxy for the bolometric luminosity or accretion 
rate of AGN (e.g., Heckman et al. 2004), since it corre- 
lates with numerous multiwavelength broad-band luminosi- 
ties (e.g., Mulchaey et al. 1994). However, since (1) the scat- 
ter of the Ln [OIII] -continuum luminosity correlation is quite 
large (see e.g., Zakamska et al. 2003, their Fig. 14), and (2) 
the Stockton & MacKenty (1987) quasar sample spans a lim- 
ited range in Ln[oiii] (3 x 10^ to 10^^ L©), the existence of the 
apparent minimum Ln[oiii] seems not to reflect a requirement 
of a minimum quasar accretion rate, instead, it may merely be 
a requirement on the availability of gaseous material in the nu- 
clear region — in order to form an EELR, a quasar may have to 
have a fair amount of gas in the nuclear NLR (given a density 
of 10^ cm-3 and [Oiii] A5007/H/3 = 10, Ln[oiii] = 6.5 x 10^^ 
erg s~^ corresponds to an ionized mass of 4.4 x 10^ Mq). In 
fact, if we replace Ln[oiii] with the absolute B magnitude of 
the quasar (Mg; compiled from Veron-Cetty & Veron 2006) 
or continuum luminosity at a certain rest-frame wavelength, 
then this correlation immediately disappears (Fig. Ic), con- 
firming that quasar accretion rate is not a dominant factor in 
the formation of an EELR. 

Combining the two correlations together, one would expect 
to achieve a high success rate in identifying EELR objects. In 
Figure Id, nuclear [O III] luminosity is plotted against radio 
spectral index for the same 47 objects for which Stockton & 
MacKenty (1987) obtained narrow-band [O III] images. If one 
pre-selects only quasars with both a^, > 0.5 and Ln[oiii] > 
6.5 X 10"^^ erg s~\ then roughly half (10/19) of the sample 
would turn out to be EELR quasars, as opposed to < 23% 
(11/47)^ 38% (10/26), or < 46% (11/24) if one performs the 
pre-selection on the basis of no a priori assumption, a^^, or 
[o III] , respectively. 

Although the identification rate of EELRs can be improved 
significantly through this kind of exercise, one inevitably en- 
counters this outstanding puzzle — why do only half of the 

^ When radio-quiet quasars are included in this kind of calculation, the de- 
rived percentage of EELR quasars should be treated as upper limits since 
Stockton & MacKenty (1987)'s survey included a disproportionally small 
number of radio-quiet quasars compared to the radio-loud ones. Note that 
there are a total of 1876 quasars known at z < 0.45 and S > -25°, among 
which less than 10% are radio-loud (Veron-Cetty & Veron 2006). 
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Fig. 1. — (a) Correlation between extended [Oiii] luminosity (Le[oiii]) and radio spectral index (au), which correlates strongly with radio morphology, (b) 
Correlation between L^iom] and nuclear [Oiii] luminosity (i-N[Oiii])- (c) ^e[OIII] versus quasar absolute B magnitude, (d) Ln[oiii] versus ajy. Open and filled 
circles represent objects with solid detections of extended [O ill] emission. Open triangles show objects where extended emission was undetected or ambiguous, 
they therefore represent upper limits of Le[oiii]- Upper limits of Ln[oiii] are drawn by arrows. Solid lines separate radio-quiet quasars (QSOs), flat-spectrum 
core-dominated quasars, and steep-spectrum lobe-dominated quasars. Dashed lines mark Lcut = 5 x 10"^^ erg s~^ which is our minimum [Oiii] luminosity for 
any extended emission to be called an EELR. So the green filled circles highlight the EELR quasars, while the open circles and triangles together are considered 
as non-EELR quasars (even though the open circles indicate detection of weak extended emission). The spectral indices of radio quiet objects (i.e., QSOs) have 
been arbitrarily set to be -0.1 for illustration purposes. Dotted lines indicate Ln[oiii] = 6.5 x 10^^ erg s~^ below which apparently no EELRs were found. 
Original data came from Stockton & MacKenty (1987) and Veron-Cetty & Veron (2006) and all except ajy have been converted to our concordance cosmology. 



Otherwise identical quasars end up having luminous EELRs 
while the other half do not? The EELR quasars and the non- 
EELR quasars not only have similar radio morphology and 
nuclear narrow-line luminosity (given that they have been pre- 
selected with both criteria), but they also look remarkably 
similar in terms of the redshift range, the broad-band quasar 
luminosity, the luminosity of the host galaxy and the central 
black hole mass. It is natural to suspect that there is at least 
one concealed parameter at work, the difference in which re- 
sults in the two distinct groups. This parameter has recently 
been identified. Quite unexpectedly, it turned out to be the 
gas metallicity of quasar broad-line regions (BLRs): quasars 
showing luminous EELRs have significantly lower metallic- 
ity BLRs (Z < 0.6 Zq) than other quasars (Z > Z©) (Fu & 
Stockton 2007a). 

2.2. The Origin of the Extended Gas 

The suggestion that the extended ionized gas is organized 
into a rotating envelope or disk (Wampler et al. 1975) was 
soon ruled out by the a counter example of 4C 37.43, where 
blue-shifted clouds were detected on both sides of the nucleus 
(Stockton 1976); instead, Stockton proposed that the ionized 
gas was ejected from the quasar itself, as strong outflows had 
been observed in a number of high-redshift quasars as blue- 



shifted absorption lines. But this scenario did not draw much 
attention. As 3C 249.1 and Ton 202 happened to be the first 
targets to be imaged in their redshifted [O III] line, and both 
showed structures similar to the galactic tidal tails in the sim- 
ulations of Toomre & Toomre (1972), Stockton & MacKenty 
(1983) argued that the extended gas is tidal debris from a 
merger of two gas-rich galaxies and the merger is also respon- 
sible for triggering the quasar activity. However, a problem 
with this scenario was noticed soon afterward. 

The density of the ionized gas was known to range from 
tens to a few hundreds cm~^, as indicated by both the 
density- sensitive [O ii] A3727 doublet profile (Stockton 1976) 
and ionization-state-sensitive [Oil] A3727/[Oiii] A5007 ra- 
tio in combination with the projected distance to the nucleus 
(Fabian et al. 1987; Crawford et al. 1988). Such a high den- 
sity implies both a small size of the ionized cloud and a high 
pressure, which is unlikely to be balanced by that of the inter- 
stellar medium of a normal galaxy. An ionized cloud of this 
density would dissipate within a sound-crossing time of ~ 10"^ 
years unless it were confined gravitationally or by the pres- 
sure of a hot surrounding medium. As there are not enough 
stars at the locations of the gaseous "tidal tails" to gravita- 
tionally confine the gas, Fabian et al. (1987) suggested that 
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the EELRs are embedded within a hot high-pressure medium. 
In this scenario, the warm ionized gas comprising the EELRs 
condenses from a hot halo sometimes seen even in poor clus- 
ters. The most direct proof for such a cooling flow would be 
a detection of diffuse X-ray-emitting gas. The X-ray lumi- 
nosity can be used to estimate the density of the hot medium 
at the locations of the EELRs (~ 10 kpc), assuming a typ- 
ical halo density profile and a typical temperature for the X- 
ray-emitting gas (e.g., = 1 keV, or T = 1 .2 x 10^ K). Unfor- 
tunately, since quasars themselves are luminous X-ray sources 
and known EELR quasars typically have a redshift of ~0.3, 
such an observation demands (1) a large- aperture X-ray tele- 
scope, (2) an arcsec-scale angular-resolution to confine most 
of the quasar emission within ~10 kpc, and (3) an accurate 
knowledge of the instrument point- spread-function to remove 
the quasar spill-over; therefore it had to wait until Chandra X- 
ray observatory was launched near the turn of the millenium. 

Contrary to the expectation, deep Chandra X-ray images of 
three EELR quasars failed to detect the hot X-ray halo from 
which the warm gas is supposed to condense, and the upper 
limits on the density of the hot gas were too low to allow suffi- 
cient cooling (Stockton et al. 2006). In fact, the pressure of the 
hot gas, as implied by the density upper limits, is comparable 
to that of the ISM in our Galaxy (nT ^ 10^ K cm"^). To re- 
iterate, without a confinement of a high-pressure surrounding 
medium, the warm ionized gas in the EELRs would dissipate 
quickly if its density is significantly greater than 1 cm~^. 

A somewhat related possibility that avoids these problems 
is that of a cold accretion flow along Mpc- scale filamentary 
structures (Keres et al. 2005). We postpone commenting on 
this option until our discussion in § 6.3. 

Following the approach pioneered by Viegas & Prieto 
(1992) and Binette et al. (1996), Robinson et al. (2000) found 
that a photoionization model of a mixed medium with both 
optically thin and thick components not only overcomes the 
outstanding difficulties in fitting certain line ratios suffered 
by previous single-phase models in the extended ionized gas 
in the radio galaxy 3C321, but it also results in an ioniz- 
ing power more consistent with far-infrared (FIR) and ra- 
dio observations of the hidden central engine where there 
is less extinction. Stockton et al. (2002) showed that the 
EELR of quasar 4C 37.43 is also better modeled as a two- 
phase medium, consisting of a matter-bounded diffuse com- 
ponent with a unity filling factor (n^2 cm~^, T ~ 15, 000 K), 
in which are embedded small, dense clouds (n ^ 500 cm~^, 
T ^ 10,000 K). Without assuming that most of the mass re- 
sides in a high-density medium, which raises the confinement 
crisis, this two-phase model can explain the high densities in- 
ferred from both the [O ll]/[0 III] ratio and the [O ll] doublet 
ratio, as the [O ll] emission is almost entirely produced by the 
dense clouds even though they account for less than 1 % of the 
total mass. 

Stockton et al. (2002) also obtained a global [O ill] velocity 
map of the 4C 37.43 EELR from spectroscopy with an image 
sheer. The velocity field appears largely disordered, consis- 
tent with the results from earlier integral field spectroscopy 
of similar objects (Durret et al. 1994; Crawford & Vanderri- 
est 1997, 2000). Some faint clouds show velocities as high as 
700 km s~\ which are unlikely to be solely gravitational in 
origin. As the embedded dense clouds would have to be con- 
stantly resupplied because of their short lifetimes, Stockton 
et al. suggested that shocks were the most likely mechanism 
for regenerating the small dense clouds. It was also known 
that two of the EELR quasars have FIR colors similar to those 



of ultraluminous infrared galaxies (ULIRGs), which are virtu- 
ally all starburst galaxies triggered by major mergers (Sanders 
& Mirabel 1996). The host galaxies of both quasars had been 
confirmed to have starburst or recent poststarburst stellar pop- 
ulations (3C 48 and Mrk 1014, Canalizo & Stockton 2000a,b). 
Combining all these pieces, Stockton et al. (2002) suggested 
that the extended emission-line gas has been expelled during 
a merger and is being shocked by a galactic superwind driven 
by a starburst triggered by the merger. 

So three decades after EELRs were first recognized, the 
origin of the ionized gas and the physical forces that con- 
trolled its distribution remained uncertain: (1) the cooling 
flow model seemed discredited following the XMM-Newton 
high-resolution spectroscopy (Peterson et al. 2003) and the 
Chandra imaging program (Stockton et al. 2006), but a new 
possibility had arisen — that the gas could have originated in 
a cold accretion flow funneled along intergalactic filaments 
(Keres et al. 2005), (2) the tidal debris model had to incorpo- 
rate a new physical mechanism (e.g., shocks from a starburst 
superwind) to explain the existence of small dense clouds em- 
bedded in the otherwise low-density photoionized medium, 
and (3) the old quasar ejection model (Wampler et al. 1975; 
Stockton 1976) had been rejuvenated because of the accumu- 
lating circumstantial evidence of shocks propagating in the 
immediate environments of EELR quasars. This evidence in- 
cludes the two-phase photoionization model of the brightest 
condensation in the EELR of 4C 37.43 (Stockton et al. 2002) 
and the discrete X-ray sources detected at distances within 40 
kpc from two EELR quasars (3C 249.1 and 4C 37.43, Stock- 
ton et al. 2006). There was thus some evidence for an out- 
flow, but it was unclear what was driving it — a starburst, or 
the quasar itself, or both? 

The ultimate goal of this and our previous papers (Fu & 
Stockton 2006; Stockton et al. 2007; Fu & Stockton 2007a,b, 
2008) has been to address these fundamental questions. The 
recent availability of integral-field units (IFUs) on large tele- 
scopes has made a dramatic improvement in our ability to 
do so. Most of the previous investigations on quasar EELRs 
have relied on either narrow-band imaging of the entire struc- 
ture or slit spectroscopy of only a small fraction of the 
whole extended emission. IFUs allow spectroscopy in a two- 
dimensional area, combining the advantages of the two tradi- 
tional observational methods. The earliest integral field spec- 
trograph used in studying quasar EELRs was TIGER (Bacon 
et al. 1995), which was a lenslet-based instrument mounted 
on the 3.6-m Canada France Hawaii Telescope (CFHT). It of- 
fered dil" xl" field of view (FOV) with a spatial sampling 
of 0.''4. Durret et al. (1994) observed three EELR quasars 
in the Stockton & MacKenty (1987) sample with TIGER (aU 
of which are included in our sample). Due to the limited 
spectral coverage of the instrument, their spectra included 
only the range from H/3 A4861 to [Olll] A5007. The main 
result of this study was that the velocity fields of EELRs 
were rather chaotic. TIGER was soon superseded by ARGUS 
(Vanderriest 1995), which was also hosted by CFHT but had 
a "lenslet+fiber" design. ARGUS covered a 8.''5 x 12.''4 
hexagonal FOV with a sampling of 0.''4. Thanks to its optical- 
fiber design, ARGUS offered a wider wavelength coverage 
than did TIGER, allowing [On] A3727 and [Oiii] A5007 to 
be observed simultaneously. Crawford & Vanderriest (1997, 
2000) performed ARGUS integral field spectroscopy of seven 
radio-loud quasars that were known to possess spatially ex- 
tended emission-line gas, among which three are in the Stock- 
ton & MacKenty (1987) sample and two are in common with 
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Fig. 2.— Narrow-band [Oiii] A5007 images of 3C48, Mrlcl014, 3C249.1, Ton616, Ton202, 4C 37.43, PKS 2251+11, and 3C79 (panels a to K). Overlaid 
are the fields of tlie CMOS IFUR observations {solid boxes) and tliose of tlie CMOS IFU2 observations (long-dashed rectangles), //^r linear-ramp-filter images 
are used if available. The radio jet direction is described by the solid lines for all of the objects except for 3C 48 and Mrk 1014, where the radio sources are less 
than 2" across. All images have been scaled to show the same physical size (150 kpc on a side) at the redshift of the object which is labeled under the scale bar. 
The images of 3C 48 and Ton 616 have been allowed to wrap around to show high surface brightness peaks. In the image of 4C 37.43, the long-dashed rectangle 
shows the mosaicked field of four IFU2 pointings, the overlapping regions of which are shown by the tick marks on the inner edge. 



our sample. Besides confirming the chaotic velocity fields, 
Crawford & Vanderriest found that all of the extended nebulae 
had very similar [O ll]/[0 III] ratios independent of location. 
These pioneering studies using IFUs on CFHT were neverthe- 
less restricted in signal-to-noise ratio (S/N) and in resolution 
because of the size of the telescope. 

3. OBSERVATIONS AND DATA REDUCTION 

The newest optical integral field spectrograph on Mauna 
Kea is the IFU (Allington- Smith et al. 2002) of the Gem- 
ini Multiobject Spectrograph (CMOS; Hook et al. 2004), 
installed on the Gemini North telescope. The IFU uses a 
"lenslet+fiber" design similar to that of ARGUS and is the 
first such instrument on an 8-m telescope. When used in 
the full-field mode, it covers a fully filled contiguous field 
of 7'' X 5'' with 1000 hexagonal lenslets 0.''2 in diameter. A 
"sky" field of 5'' x 3. ''5 V from the science field is observed 
simultaneously by another 500 lenslets for background sub- 
traction. Optical fibers behind the lenslet arrays reformat the 
focal plane into two parallel pseudo-slits that pass the light 
into the rest of the spectrograph. To avoid spectral overlap- 
ping between the two pseudo-slits, broad-band filters are of- 
ten used to limit the spectral range. One can increase the 
spectral coverage to that of long-slit spectroscopy by sacri- 
ficing half of the FOV (i.e., masking off one of the slits). This 
instrument provides sufficient spatial and spectral resolution 
to allow constructing unprecedentedly high-resolution veloc- 
ity fields and resolving the density- sensitive [O ll] doublet at 
redshifts beyond ~0.3. In addition, the wide wavelength cov- 
erage in combination with the large aperture of Gemini permit 
detailed spatially resolved modeling of emission-line spec- 
tra accounting for weak but important lines such as [Ne V] 
A3426, He ii A4686 and [O ill] A4363 that were unexplored 
by previous IFU observations. 

The narrow-band [Olll] imaging survey of Stockton & 



MacKenty (1987) pubHshed [Olll] A5007 luminosities and 
morphologies for all cases of detected extended emission 
around 47 quasars with z < 0.45, S > -25°, and My < 
-23.4^. We have performed GMOS integral field spec- 
troscopy of 7 of the 11 quasars that show luminous struc- 
tures of extended emission with Le[oiii] > 5 x 10"^^ erg s~^ — 
3C 249.1, 4C 37.43, 3C48, Mrk 1014, Ton 616, Ton 202, and 
PKS 2251+11. We also have obtained comparable observa- 
tions of a radio galaxy, 3C 79, which has a similarly luminous 
EELR. The properties of these objects have been listed in Ta- 
ble 1. The literature redshifts (z') and absolute B magnitudes 
(Mb) in this table are compiled from Veron-Cetty & Veron 
(2006). The more accurate redshifts (z) are measured from the 
[Oiii] AA4959,5007 lines in nuclear spectra extracted from 
our GMOS/IFU data. The radial velocity maps presented 
in § 4 have been calculated relative to these assumed sys- 
temic redshifts using the relativistic Doppler formula (Stock- 
ton et al. 2007). In addition to the rest-frame correction, the 
velocity dispersion (a = FWHM/2.355) maps have also been 
corrected for the instrumental resolution. The typical spec- 
tral resolutions of the B600 and the R400 gratings are re- 
spectively about 2.1 A and 2.7 A FWHM, as measured from 
night sky lines. The last four columns give the continuum lu- 
minosity at H/3, nuclear [O ill] A5007 luminosity, extended 
[Oiii] luminosity within an annulus of inner radius 11.2 kpc 
and outer radius 44.9 kpc centered on the nucleus, and ra- 
dio spectral index between 400 and 2700 MHz. Results for 
the quasars 3C 249.1 and 4C 37.43, and for the radio galaxy 
3C79, have been pubHshed in Fu & Stockton (2006, 2007b, 
2008). In this paper we present results for the last five ob- 

^ This absolute V magnitude cut was calculated in their original cosmology 
(Hq = 75 km s~^ Mpc~^ , Q,n = 0, and Oa = 0) without /^-correction. In our 
chosen cosmology and accounting for /^-correction, all the 47 objects have 
Mb < -22.4 (see Fig. Ic). 
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TABLE 1 

Targets of Gemini GMOS Integral Field Spectrocopy 



Name 


Designation 


z' 


z 


Mb 


^A4861 


^N[OIII] 


^E[OIII] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C 48 


0137+3309 


0.367 


0.36926 


-24.59 


41.29 


43.22 


42.75 


0.82 


Mrk 1014 


0159+0023 


0.163 


0.16323 


-23.14 


41.08 


42.88 


41.81 


RQ 


3C 249.1 


1104+7658 


0.313 


0.31172 


-25.18 


41.71 


43.43 


42.82 


0.84 


Ton 616 


1225+2458 


0.268 


0.26808 


-23.17 


40.84 


42.94 


42.59 


0.80 


Ton 202 


1427+2632 


0.366 


0.36367 


-25.39 


41.37 


43.12 


42.32 


0.62 


4C 37.43 


1514+3650 


0.370 


0.37120 


-24.09 


41.45 


43.28 


42.97 


0.80 


PKS 2251+11 


2254+1136 


0.325 


0.32538 


-24.88 


41.36 


42.94 


42.47 


0.73 


3C 79 


0310+1705 


0.256 


0.25632 






42.21 


42.89 


0.87 



Note. — Col. (1): Common name. Col. (2): J2000.0 designation. Col. (3): Redshift from 
literature. Col. (4): Redshift measured from our GMOS/IFU data from the [Oiii] AA4959,5007 
lines. Col. (5): Absolute B magnitude (/:-correction appHed). Col. (6): Continuum luminosity at Y{(3 
in logarithmic (erg s~^ A~^). Col. (7): Nuclear [Oiii] A5007 luminosigy in logarithmic (erg s~^). 
Col. (8): Extended [Oiii] A5007 luminosity in logarithmic (erg s~^). Col. (9): Radio spectral 
index (^ oc between 400 and 2700 MHz. Data in Cols. (3) & (5) are taken from Veron- 

Cetty & Veron (2006). Cols. (6-9) are based on Stockton & MacKenty (1987). All luminosities 
and magnitudes in this table have been converted to our concordance cosmology. 3C 48 is a CSS 
radio-loud quasar with a single-sided jet within r^l" to the north (Feng et al. 2005); Mrk 1014 is a 
radio-quiet quasar but shows two radio knots r^\" to either side of the nucleus (Leipski et al. 2006; 
Bennert et al. 2002) and has a steep radio spectrum {ai, ~ 0.9); 3C 79 is an FR-II radio galaxy. The 
rest of the objects are all FR-II radio-loud quasars. 



TABLE 2 

Gemini GMOS Integral Field Spectroscopy Observing Log 



Slit 


Grating 


CW (A) 


Lines 


RA. Field Exposure (s) 


Airmass 


Seeing 


UT Date 


(1) 


(2) 


(3) 


(4) 


(5) (6) (7) 


(8) 


(9) 


(10) 


3C 48 (z = 0.369) 


IFUR 


B600 


5750 


[NevHOiii] 


30° on 1800x5 


1.201 


W.'A5 


20030923 


IFUR 


B600 


5680 




90° off 2700x3 


1.468 


O'.'ll 


20070819 


IFU2 


R400 


9000 


[OiHSii] 


90° on 2400x1 


1.094 




20070905 


Mrk 1014 (z = 0.163) 


IFUR 


R400 


6250 


H^-[S II] 


90° on 2400x1 


1.114 


o:^54 


20070813 










off 2400x1 


1.070 






IFUR 


B600 


4890 


[Nev]-[Oiii] 


90° on 2400x1 


1.245 


0fM9 


20070923 










off 2400x1 


1.403 


0'.'13 


20070924 










on 2400x1 


1.430 




20071007 










off 2400x1 


1.618 


Of'50 


20071012 


3C 249.1 (z = 0.312) 


IFUR 


B600 


5800 


[NevHOiii] 


130° on 2825x2 


1.841 


i;^15 


20050407 










on 2825x1 


1.838 


lf'03 


20050408 


IFUR 


R400 


7170 


H/3-[S II] 


130° on 1575x1 


1.840 


I'.' 03 


20050408 


Ton 616 (z = 0.268) 


IFUR 


R400 


6570 


H^-[S II] 


60° off 2400x3 


1.207 


0'.'51 


20070612 


IFUR 


R600 


7250 


H^-[S II] 


60° on 1800x3 


1.757 


O'.'IA 


20070717 


IFU2 


B600 


6332 


H^-[0 III] 


0° on 2400x1 


1.230 


0'.'90 


20070703 


Ton 202 (z = 0.364) 


IFUR 


B600 


5680 


[Nev]-[Oiii] 


340° on 2100x3 


1.108 


O'.'^l 


20070414 










off 2100x3 


1.197 






IFU2 


R600 


9107 


[Oi]-[Sii] 


340° on 1800x2 


1.069 


Of'58 


20070709 


4C 37.43 (z = 0.371) 


IFUR 


B600 


6412 


[Nev]-[Oiii] 


0° on 2400x5 


1.071 


o:^59 


20060524 


IFU2 


R831 


6585 


H^-[0 III] 


0° on 720x12 


1.059 


0i'A9 


20060523 


PKS 2251+11 (z = 0.325) 


IFUR 


B600 


5510 


[Nev]-[Oiii] 


217° on 2700x3 


1.064 


0'.'39 


20070811 










off 2700x2 


1.031 














off 2700x1 


1.059 


0'.'33 


20070812 


IFU2 


R400 


8810 


[Oi]-[Sii] 


217° on 2400x1 


1.224 


0'.'39 


20070905 


3C 79 (z = 0.256) 


IFUR 


B600 


5500 


[Nev]-[Oiii] 


160° on 2880x3 


1.253 


O'.'Sl 


20061221 



Note. — Col. (1): IFU configuration (IFUR: half-field one-slit; 1FU2: full-field two-slit). Col. (2): Grating. Col. (3): Central 
wavelength. Col. (4): Emission lines enclosed (Abbreviations refer to [Ne v] A3426, [O III] A5007, [O i] A6300, and [S ii] A6731). 
Col. (5): Position angle. Col. (6): Pointing position ("on-nucleus" or "off-nucleus", refer to Fig. 2), if there are multiple pointings. 
Col. (7): Exposure time. Col. (8): Mean airmass during the exposure. Col. (9): Seeing measured from acquisition images. Col. 
(10): Observation date (yyyymmdd). "-" means "the same as above". 
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jects. The observations were performed in queue mode in 
2007 between April and October. Various slit/grating/central 
wavelength combinations were used to avoid chip gaps for 
objects with different redshifts. Generally the single- slit half- 
field mode (IFUR) was used to cover the wider wavelength 
range from [Ne V] A3426 to [Oiii] A5007 (blue spectra), and 
the two-slit full-field mode (IFU2) was used to cover the nar- 
rower wavelength range between [O i] A6300 and [S ii] A673 1 
(red spectra). A detailed observing log can be found in Ta- 
ble 2. Figure 2 shows the areas covered by the 3f'5x5'' FOV 
of the IFUR mode and the l"x5" FOV of the IFU2 mode. 
Most of the observations were performed under excellent see- 
ing (FWHM < O'.'G). Spectrophotometric standard stars were 
observed with the same instrument settings as that of the sci- 
ence observations but often on different nights for relative flux 
calibration. 

The data were reduced using the Gemini IRAF package 
(ver. 1.8), following the same procedure as described by Fu 
& Stockton (2007b). After correction for differential at- 
mospheric refraction, the quasars were removed from the 
data cubes using the same Richardson-Lucy method as in 
Fu & Stockton (2007b) with the STSDAS procedure cplucy 
and point-spread functions (PSFs) generated in emission- 
line-free regions in the same data cube. 

To increase the S/N, we extracted spectra of distinct clouds 
with circular apertures of a radius between 0.''5 to 0.''7, de- 
pending on the seeing. We have concentrated on clouds at 
least 1 .''5 away from the nucleus to avoid high levels of con- 
tamination from quasar light to the emission lines. The line- 
of- sight Galactic reddening was corrected using a standard 
Galactic reddening law (Cardelli et al. 1989) with extinction 
values from Schlegel et al. (1998). Intrinsic reddening due 
to dust in front of or associated with the cloud was deter- 
mined from the Ha/H/3 or H7/H/3 ratio, assuming theoretical 
ratios of Ha/H/3 = 3.1 and H7/H/5 = 0.468 (Osterbrock 1989). 
Whenever Ha and H/3 were observed simultaneously, the 
Ha/H/3 ratio was used to determine the intrinsic reddening; 
and in other cases, the H7/H^ ratio was used. The emission- 
line fluxes were dereddened using the same reddening law and 
are tabulated in Table 3. When there are overlaps between 
blue and red spectra, which is the case for Mrk 1014, the spec- 
tra were normalized at the overlapping region according to the 
[Oiii] AA4959,5007 flux and combined before dereddening. 
For 3C48, Ton 202 and PKS 2251+11, as there are no over- 
laps, the dereddened line intensities from the red spectra were 
divided by a constant so that Ha/H/? = 3.1. In contrast to the 
case of the intrinsically weaker extended emission in Seyfert 
galaxies, we do not here have to deal with the issue of cor- 
recting for stellar B aimer absorption in measuring the B aimer 
emission lines to determine the internal extinction. Most of 
the EELRs we are dealing with are projected onto regions 
with no evidence for stellar contamination detectible in our 
observations. Even for those in the inner regions, where we 
can detect the underlying host galaxy, the emission so strongly 
dominates the stellar continuum that absorption corrections 
are negligible. 

4. RESULTS 
4.1. 3C48 

3C48 (z = 0.369, I" = 5.11 kpc) was not only the first 
quasar to be identified (Sandage & Matthews 1961; Green- 
stein & Matthews 1963), but was also the first quasar to be 
associated with a luminous EELR (Wampler et al. 1975). It 



is a compact steep- spectrum (CSS) radio-loud quasar with a 
one-sided jet appearing disrupted within ~0.''7 (3.6 kpc) to 
the north (Wilkinson et al. 1991; Feng et al. 2005). 

The host galaxy of 3C 48 is unusually luminous among low- 
redshift quasars (Kristian 1973). Strong Balmer absorption 
lines were detected long ago in the host galaxy, indicating the 
existence of a relatively young stellar population (Boroson 
& Oke 1982). With deep spatially-resolved long-slit spec- 
troscopy and a stellar population decomposition technique, 
Canalizo & Stockton (2000a) determined the ages for the 
most recent major star formation in 32 distinct regions within 
50 kpc of the nucleus, and they detected ongoing star forma- 
tion and post-starburst populations less than 115 Myr old in 
26 of the 32 regions. Consistent with the spectroscopy results, 
the IRAS FIR luminosity and colors of 3C 48 look quite simi- 
lar to those of ULIRGs (Stockton & Ridgway 1991); virtually 
all of the latter appear to be merger- or interaction-induced 
starbursts (Sanders & Mirabel 1996). The optical morphol- 
ogy of the host galaxy also indicates an ongoing merger. 
A tidal tail comprising an old stellar population is seen ex- 
tending to the northwest (Canalizo & Stockton 2000a), and 
there is a bright continuum peak northeast of the nucleus 
(3C48A; Stockton & Ridgway 1991) comprising stars with 
a luminosity-weighted age of 140 Myr. 3C48A seems to be 
the distorted nucleus of the merging partner, and its stellar 
age may indicate star formation triggered during the previous 
close encounter (Stockton et al. 2007). 

3C 48 was first imaged in [O III] through a narrow-band fil- 
ter by Stockton & MacKenty (1987). Like the radio jet, the 
bulk of the emission-line gas is distributed to the north of the 
quasar nucleus (see Fig. 2a for a HST [O III] linear-ramp fil- 
ter [LRF] image, Kirhakos et al. 1999). Wampler et al. (1975) 
noticed that the [O III] A5007 line profile is unusually broad — 
26 A FWHM in rest-frame (i.e., - 1560 km s"^). They also 
found strong Fell emission between H7 and ¥1(3. Subse- 
quent higher resolution spectra showed that the [O III] line is 
double-peaked (Chatzichristou et al. 1999; Canalizo & Stock- 
ton 2000a): in addition to the usual nuclear NLR at the same 
redshift as the broad HP line, there is an anomalous high- 
velocity narrow-line component 0.''25 north of the nucleus, 
precisely in the initial direction of the radio jet, but within the 
region where the jet appears relatively unperturbed (Stock- 
ton et al. 2007). Chatzichristou et al. (1999) also presented 
an [O III] A5007 velocity map covering a large fraction of the 
EELR, and they found no clear trends in the velocity field ex- 
cept for an apparent rotation within V to the nucleus. 

In addition to the IFUR data cube that was used to study 
the high- velocity anomalous cloud and the 3C 48A (both are 
within 1'' from the nucleus), we obtained another IFUR data 
cube covering the EELR clouds further north and an IFU2 
data cube covering the whole EELR (Fig. 2a and Table 2). 
As the IFU2 observation has the most complete spatial cov- 
erage and an excellent resolution, ideally we should construct 
a global velocity field from the IFU2 data cube alone using 
the Ha line, which is the strongest line covered by the cube. 
However, the quasar was placed very close to the edge of the 
IFU2 FOV, making it difficult to remove the quasar scattered 
light. Therefore, we built a hybrid data cube by combining 
the [O III] A5007 region in the on-nucleus IFUR cube with the 
Ha region in the IFU2 cube. First, the quasar was carefully 
removed using the Lucy-Richardson technique from the IFUR 
cube. Second, we shifted and interpolated the IFU2 Ha cube 
to accurately match the velocities of the IFUR [O III] cube. 
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Fig. 3. — Channel map and velocity field of the EELR of 3C48. The plus signs indicate the position of the quasar before it was removed from the data cube. 
North and east are marked by the directional arrow. Scale bars of 5 kpc in length at the quasar redshift are shown. The channel maps are shown above the velocity 
field which covers the same area. The central velocities, in km s~^ relative to that of the nuclear NLR (refer to Table 1) are indicated in the top right corner of 
each panel. The tick marks on the axes are spaced by \" . Regions where full spectra have been extracted are labeled. The three columns of the velocity field 
are line intensity (*S), radial velocity (V), and velocity dispersion (cr) maps. The anomalous cloud 0f^25 north of the nucleus has V ~ -500 km s~^ and a ~ 
390 km s~^ and therefore exceeds the displayed range of the maps. Contours are from the intensity map, where the green arrow emanating from the nucleus 
delineates the compact radio jet of 3C 48 (the arrow ends at the component N2a in Feng et al. 2005). The data outside of the dash-dot rectangle are from the Ha 
line in the IFU2 data cube, which has been scaled and interpolated to match both the flux and the velocities of the [O ill] A5007 line in the IFUR cube inside the 
rectangle. 



then we scaled the former by a constant to match the fluxes 
of the clouds 1" north of the quasar, which were covered by 
both cubes, to those of the IFUR cube. Finally, we rotated the 
IFUR cube and inserted it into the IFU2 cube in the overlap- 
ping area of the two. Figure 3 shows the channel maps and the 
velocity field from this hybrid data cube. Our velocity maps 
are in generally good agreement with the previous integral 
field spectroscopy of Chatzichristou et al. (1999)^, but they 
offer a much finer spatial and velocity resolution. The highest 
relative velocities were observed in the anomalous cloud 0.''25 
to the north and the group of 3 clouds ~ 1" also to the north. 
The three-cloud group is also seen in the original IFU2 data 
cube in Ha at the same relative velocities but with much lower 
S/N. The high- velocity cloud near the nucleus is also the only 
region where high velocity dispersions (a ^ 100 km s~^) are 
seen. 

We extracted a nuclear spectrum from the on-nucleus IFUR 
cube with an 0. ''5 -radius aperture centered on the quasar, from 
which we measured the redshifts of the high- velocity cloud 
and the NLR at systemic velocity. We found that (1) the 

^ Note that their relative velocities were calculated with respect to an arbi- 
trary systemic velocity at z = 0.368. 



high- velocity and systemic components were at redshifts z - 
0.36674±0.00006 and z = 0.36926±0.00002; therefore there 
is a velocity difference of -501 ±16 km s~\ and (2) their ve- 
locity dispersions (a) are 392±21 and 148±5 km s~^ The 
above nominal values and their errors were respectively the 
mean values and half of the differences measured from the 
[Oiii] A5007 line and the [Oiii] A4959 line. There are slight 
differences between our results and those of Stockton et al. 
(2007), because we corrected for the atmospheric 5-band ab- 
sorption with a standard star spectrum taken on the same 
night. 

The [O II] AA3726, 3729 doublet is wefl-resolved in the two 
brightest clouds (c and d). Since their blue spectra come from 
the off-nucleus IFUR data cube which was taken during a 
period of poorer seeing (~ 0.''7), the two clouds were not 
resolved as clearly as in the IFU2 data cube. So, we com- 
bined their spectra and measured that [Oil] A3726/A3729 
= 0.946±0.036 and [Oiii] (A4959 + A5007)/A4363 = 51±4, 
which implies an electron density (Ne) of 360±50 cm~^ and 
an electron temperature (7;) of 17400±700 K (STSDAS task 
temden). 
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Fig. 4. — Channel map and velocity field of the EELR of Mrk 1014. Keys are the same as in Fig. 3. Note that (1) there is a large gap in velocity between the 
first row and the second row of the channel map, and (2) the emission-line cloud appeared 1" to the west of the nucleus in panels at -1177 and -1105 km s~^ 
(This cloud is not shown in the velocity field, as it spatially overlaps with a more luminous cloud at ~150 km s~^). The dash-dot circles show the two off-nuclear 
radio components — they approximate the 0.06 mJy/beam contour in a VLA 8.4 GHz map (Leipski et al. 2006). 



4.2. Mrkl014(PG0157+001) 

Mrk 1014 {z = 0.163, V' = 2.80 kpc) was discovered as 
an ultraviolet excess object in the Byurakan objective prism 
survey (Markarian et al. 1977). Subsequent spectroscopy 
observed a broad 11(3 line, confirming its nature as a low- 
redshift radio-quiet quasar (Afanasev et al. 1980) with an 
radio-loudness (R) parameter of 2.1 (Kellermann et al. 1989). 
The low-luminosity radio source was resolved into a triple 
source at a resolution of Of' 3 6 by the Very Large Array 
(VLA) — two radio knots within l.'^l (3.1 kpc) on either side 
of a central component that corresponds to the nucleus (Leip- 
ski et al. 2006; Miller et al. 1993; Kukula et al. 1998). The 
host galaxy of Mrk 1014 is very similar to that of 3C 48. The 
FIR luminosity and colors qualify Mrk 1014 as a ULIRG. It 
has a prominent spiral-like tidal tail extending to the north- 
east, the spectra of which showed stellar absorption features 
of a "K+A" character (MacKenty & Stockton 1984). Strong 
Fe II emission is also present in the nuclear spectrum. Post- 
starburst populations of ages between 1 80 and 290 Myr were 
detected in 5 distinct regions in the host galaxy, which, to- 
gether with its morphology, suggest that Mrk 1014 is a result 
of a merger of two disk galaxies of comparable size (Canalizo 
& Stockton 2000b). 

Mrk 1014 is the only radio-quiet quasar with a luminous 
EELR in the Stockton & MacKenty (1987) sample (Fig. la). 
The morphology of the EELR, although on a much larger 
scale, resembles that of the radio emission (Fig. 2b). There 
is evidence for strong jet-cloud interactions within the extent 
of the radio components (< l.'^l or 3.1 kpc). By subtracting 
an /f^T narrow-band line-free image off an [O III] A5007 LRF 
image, Bennert et al. (2002) identified a filamentary emission- 
line structure 0.''9 to the west of the nucleus, which appeared 
coincident with the western radio knot 1 1 from the central 
radio component at a RA. of -82.2°. Near the same region, 
Canalizo & Stockton (2000b) detected a weak high- velocity 
cloud at -1350 km s~^ as well as a strong emission-line cloud 



at -200 km s~^ (both velocities were calculated with respect to 
the underlying stellar features). This high- velocity cloud was 
confirmed by Leipski & Bennert (2006), who also noticed that 
there was a broad (a ~ 590 km s~^) blueshifted component 
in addition to the classical NLR in the nuclear [O III] lines, 
which are reminiscent of the double-peaked [Olll] lines in 
3C48. 

When removing the quasar nucleus from the data cube with 
plucy, Mrk 1014 had to be dealt with differently from the 
rest of the sample because the quasar is not the only non- 
negligable continuum source near the nucleus. Bennert et al. 
(2002) pointed out that the HST line-free continuum image of 
Mrk 1014 was clearly extended compared with that of a star. 
We also observed a significant decrease in the FWHM of the 
central source at the wavelengths of the broad H/? and H7 
lines in the data cube, implying the presence of a continuum 
source that has a spectrum different from that of the quasar. 
So, we decided to use the on-nucleus IFUR data cube of 3C 48 
to create the PSF because it was observed with a very similar 
instrument setting and under similar seeing. Due to the well 
known wavelength dependence of the spatial resolution, we 
had interpolated the 3C 48 data cube to the same wavelength 
grid as that of the Mrk 1014 cube before generating the PSF. 

Figure 4 shows the channel maps and the velocity field af- 
ter removing the quasar nucleus. A peanut- shaped continuum 
source is evident at all velocities. It is unclear what the na- 
ture of the continuum source is from the current data, but it is 
possibly synchrotron radiation from the radio jets connecting 
the radio core and the knots. The high- velocity cloud I" to the 
west is seen in the channel maps centered at -1 177 and -1015 
km s~^and appears offset from the radio knot to the southeast. 
At almost the same location, another emission-line cloud is 
seen at lower velocities which appears to have a much higher 
velocity dispersion and to be offset from the radio knot to the 
northeast (see the panel at +277 km s~^). Although rather dis- 
ordered, the velocity field shows that most of the line-emitting 
clouds have velocities within 200 km s~^ from the systemic 
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Fig. 5. — Channel map and velocity field of the EELR of Ton 616: central region. Keys are the same as in Fig. 3. To separate multiple clouds that are present 
along the same line of sight, the velocity field is shown in two velocity bands (-170 to 170 km s~^ and 105 to 470 km s~^). 



velocity. Another similarity to the EELR of 3C 48 is that the 
only region of high velocity dispersions (a ^ 100 km s~^) is 
the area occupied by the radio structure. 

Although the [O ll] doublet is unresolved as a result of the 
low redshift of Mrk 1014, the [S ii] AA6716, 6731 doublet al- 
lows accurate measurements to be made. This is possible be- 
cause at z < 0.26 the [S ll] lines remain in the /'-band where 
CCD fringing is much less severe than in the z'-band. Unfor- 
tunately, the latter is the case for the other four quasars. The 
ratio of the [S ll] doublet can be used to infer electron densi- 
ties (Osterbrock 1989). We measure [Sii] A6716/A6731 ra- 
tios of 1.30±0.03, 1.29±0.03, 1.30±0.03, and 1.52±0.09 in 
regions a, b, c, and d, respectively. The results for a, b, and c 
imply a density of A^, = 121±30 cm"^ if 7; = 10^ K, or A^, = 
1 17±40 cm-^ if = 1.5 X 10^ K. The [S ii] ratio of region d 
is at the low density limit: Ne < 80 cm~^ . 

The atmospheric A-band absorption affects the intensities 
of the [Nil] A6548, Ha and [Nil] A6583 lines; therefore their 
intensities had to be corrected. The empirical absorption 
curve that we used to correct for the A-band absorption was 
derived from an observation of an AOV star (HIP 13917) taken 
at a similar airmass immediately after the science exposures. 



Ton 616 (z = 0.268, V = 4.11 kpc) was first cataloged as a 
Tonantzintla blue star (Iriarte & Chavira 1957). The bright ra- 
dio source linked with Ton 616 was found in the Fourth Cam- 
bridge radio survey (4C; Pilkington & Scott 1965). Subse- 
quent spectroscopy identified Ton 616 as a quasar with a red- 
shift of 0.268 (Lynds & Wills 1968). The radio emission has 
a classical FR II type morphology, showing hot spots ~35'' 
(144 kpc) to the northeast and southwest of a faint core at a 
PA. of 30° (Gower & Hutchings 1984). 

The host galaxy of Ton 616 is an elliptical galaxy with an 
effective radius of ri/2 — 20 kpc^, as measured in a ground- 
based /T-band image (Dunlop et al. 1993; Taylor et al. 1996). 
No signs of disturbance have been reported in the literature. 
Unlike those of 3C 48 and Mrk 1014, the nuclear spectrum of 
Ton 616 does not show strong Fell emission (e.g., Baldwin 
1975; Steiner 1981). 

The EELR of Ton 616 consists of (1) an inversion- 
symmetric distribution of gas in the vicinity of the quasar 
nucleus, (2) a detached arch 7. ''5 to the north, (3) and some 
fainter clouds scattered in the south and southwest (Stock- 
ton & MacKenty 1987). Our integral field spectroscopy com- 
pletely covered the first two components (Fig. 2d). The veloc- 
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Fig. 6. — Channel map and velocity field of the EELR of Ton 616: the northern arch. 



ity field of the central region and the northern arch are shown 
in Fig. 5 and 6, respectively. Partly due to the mediocre see- 
ing (0.''9) of the IFU2 observations, line blending (i.e., multi- 
ple kinematic components present in the [O III] A5007 line of 
a single spectrum) occurs quite often. Thus, we display the 
velocity field in two velocity bands. The three clouds com- 
prising the northern arch are essentially at the same velocity 
(^ -160 km s~^). Our results are in good agreement with 
the velocity map of Durret et al. (1994) from TIGER inte- 
gral field spectroscopy and the rather coarse map of Boisson 
et al. (1994) from long-slit spectroscopy. The much higher 
spectral resolution of GMOS/IFU in comparison with TIGER 
allowed us to cleanly disjoin spatially overlapping but kine- 
matically independent components in the two regions where 
Durret et al. noticed line widths significantly broader than the 
instrumental resolution. 

Since only the southwest part of the central IFU2 field was 
observed with a wide spectral coverage setting (i.e., the on- 
nucleus IFUR observations), we could only extract spectra 
from two clouds (a and b) for the emission-line analysis. Un- 
fortunately, the observations meant to cover the blue spec- 
tra were never executed — we thus do not have intensity mea- 
surements for emission lines shortward of He II A4686. The 
data cube of the northern arch in fact covers a wavelength 
range from [O ii] A3727 to [S ii] A6731. But as a result of the 
low efficiency « 40 %) of the R400/G5305 grating at A < 
5500 A, all of the emission-lines at wavelengths shorter than 
He II A4686 are too noisy to be measured. 

Only in region e is the [S ll] doublet sufficiently strong rela- 
tive to the noise to allow a meaningful measurement of its ra- 
tio: [Sii]A6716/A6731 = 1.33±0.14, implying A^, = 90±160 
(80±200) cm-3 if Te = 10^ (1.5 x 10^) K. 

4.4. Ton 202 (B2 1425+267) 

Ton 202 (z = 0.364, \" = 5.07 kpc) was also one of the 
Tonantzintia blue stars (Iriarte & Chavira 1957). The ra- 
dio source was first listed in the B2 catalogue (Colla et al. 
1972). Ton 202 was first suspected to be a quasar by Green- 
stein (1966) based on an ambiguous detection of several broad 
emission lines, and this identification was subsequently con- 
firmed and the object was assigned a redshift of 0.366 as mea- 
sured from the broad ¥1(3 line (Weedman 1976). The radio 
source has an FR II morphology with one radio hot spot 140^' 



(709 kpc) to the northeast at a R A. of 52° and another one 83'' 
(420 kpc) to the southwest at a RA. of -120° (Rogora et al. 
1986). 

By modeling an HST WFPC2 broad-band images of 
Ton 202, Kirhakos et al. (1999) classified the host galaxy as 
an elliptical and measured an effective radius of r 1/2 — 12 
kpc. An arc-like feature to the west of the nucleus was seen 
in the HST F555W image (Kirhakos et al. 1999) and a CFHT 
//-band adaptive optics image (Marquez et al. 2001), suggest- 
ing mild galactic interactions. This feature was also detected 
by Stockton & MacKenty (1983) in their line-free continuum 
image, but the feature is not seen in the HST F814W image. 
It is clear from these images that the host galaxy of Ton 202 
is not experiencing as dramatic interactions as those of 3C 48 
and Mrkl014. Like that of Ton 616, the nuclear spectrum 
of Ton 202 does not show strong Fell emission (e.g., Steiner 
1981;Boroson etal. 1985). 

The EELR of Ton 202 was discovered by Stockton & 
MacKenty (1983) in their narrow-band [Olll] imaging sur- 
vey. The nebulosity shows two impressive spiral-like arms 
extending over 40 kpc (Fig. 2e). Long-slit spectroscopy in the 
[O III] showed that the brighter eastern arm was approaching 
at a velocity of -170 km s~^ while the fainter western arm is 
receding at +190 km s~^ with respect to the systemic velocity 
(Boisson et al. 1994). We present in Fig. 7 the [O ill] channel 
maps and the velocity field of the brighter eastern arm from 
the mosaicked IFUR data cube. Similar trends were seen in 
Ha, from the IFU2 data cube (which only covers the Ha - 
[S II] region), although the data are much noisier. We mea- 
sured a velocity range of -250 to -100 km s~^ in the bright- 
est condensation a, consistent with the velocities reported by 
Boisson et al. (1994). The entire velocity field is definitely 
more complex than pure rotation. 

Because of the poor S/N at both the [O ii] and the [S ii] 
doublets, we were unable to measure the electron density in 
this EELR. 

4.5. PKS2251+11 

PKS 2251+11 (z = 0.325, \" = 4.70 kpc) was discovered 
in the Parkes radio sky survey (Day et al. 1966). Soon the 
optical counterpart was identified (Bolton & Ekers 1966), and 
the object was confirmed to be a quasar from its ultraviolet 
excess (Kinman et al. 1967) and a redshift was determined 
(Kinman & Burbidge 1967). The extended radio emission 
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Fig. 7. — Channel map and velocity field of the EELR of Ton 202. Keys are the same as in Fig. 3. 



'PKs'225i + 1l' -479 
[0 III] km S-' 

- ^ s/- 


-363 


" -247" 


" -132" 






+ 214 

* : 






+561" 

^ - 



S (10"^^ erg s"^ cm"2 orcsec"^) V (km s"^) 

12 3 4 5 6 -300 -200 -100 100 200 300 




Fig. 8. — Channel map and velocity field of the EELR of PKS 2251+11. Keys are the same as in Fig. 3. The white dash dot ellipses in some panels describe 
the radio hot spots — they approximate the 12.7 mJy/beam contour in a VLA 4.86 GHz radio map (Price et al. 1993). 



is an FR II source with a P. A. of 138° and an total extent 
of only 11'' (52 kpc), and the southeast lobe appears much 
closer to the core than the northwestern lobe (Price et al. 1993; 
Miller et al. 1993). Optical broad-band imaging revealed that 
PKS 2251+11 is hosted by an elliptical galaxy with no signs 
of disturbance (Hutchings & Neff 1992). A moderate level of 
Fe II emission is present in the nuclear spectrum of the quasar 
(e.g., Hutchings & Crampton 1990). 



Stockton & MacKenty (1987) noticed that there was a clear 
correspondence between the southeast part of the EELR and 
the radio structure based on a lower resolution VLA map 
(Hintzen et al. 1983). One of the two bright emission-line 
clouds ~ A" southeast of the quasar appears to be coinci- 
dent with a hot spot in the southeast radio lobe (Fig. 2g). 
Hutchings & Crampton (1990) caught a glimpse of the chaotic 
velocity structure of the EELR by taking off-nucleus spec- 



Quasar Extended Emission-Line Regions 



13 




-0.5 0.0 0.5 

Logio([0 II] 3727/[0 III] 5007) 




0.5 1.0 1.5 

Log,o([0 II] 3727/[Ne III] 3869) 



> 0.5 




-2.0 



- 1 .5 - 1 .0 

Logio([0 I] 6300/Ha) 




-0.5 

Logio([S II] 6717,31/Ha) 



Fig. 9. — Line ratios of quasar EELRs indicate photoionization by an AGN-type ionizing spectrum. Log-scaled density distributions of SDSS DR4 
emission-line galaxies are displayed in line-ratio diagrams: [O ill] X5001/Uf3 vs. (a) [O ii] AA3726, 3729/[0 ill] A5007, (b) [O ii] AA3726, 3729/[Ne ill] A3869, 
(c) [O I] A6300/Ha, and (d) [S ii] AA6717, 6731/Ha. Measurements from various clouds in the EELRs of 3C 48, Mrk 1014, Ton 616, Ton 202, and PKS 2251-^1 1 
are shown as color symbols with error bars. The symbol colors and shapes follow the coding system in the legend (e.g., 3C 48-^ is indicated by the purple circle). 
Seyfert2s, LINERs and star-forming galaxies are well separated, by definition, in the last diagram. However, they appear quite blended in the first three diagrams. 
Therefore, Seyfert2s are shown as blue contours, LINERs as red contours, and star-forming galaxies (including star-forming — AGN composite galaxies) as the 
background image. Arrows are reddening vectors. 3C48-fl is the only cloud that has a significant star-forming component (see also Fig. lla). 



tra along three slit positions. IMore detailed velocity maps 
became available as the object was investigated by integral 
field spectroscopy (Durret et al. 1994; Crawford & Vanderri- 
est 2000). A bipolar velocity structure was explicitly pointed 
out by Crawford & Vanderriest (2000) — redshifted clouds are 
distributed in an envelope in the vicinity of the nucleus, while 
more extended clouds are mostly blueshifted. It was also re- 
alized that a physical connection between the southeast cloud 
and the radio jet based on their spatial correspondence might 
be problematic, since the line width of the emission-line cloud 
is rather narrow (a < 170 km s~^). 

Our GMOS data at 0'/4 resolution are presented as chan- 
nel maps and two-dimensional velocity fields in Fig. 8. The 
results are in general agreement with previous observations 
(Hutchings & Crampton 1990; Durret et al. 1994; Crawford & 
Vanderriest 2000). In a couple of the channel maps, we have 
overlaid the radio structure on top of the optical emission-line 
clouds. This comparison shows that the southeast cloud (a) 
that is alleged to coincide with the radio peak is actually a bit 
off-centered with respect to the hot spot. In agreement with 
Crawford & Vanderriest (2000), no increase in velocity dis- 
persion is seen in cloud a relative to the rest of the EELR, 
indicating that the connection between the radio plasma and 
the warm ionized gas is not a physical one (i.e., it is most 
likely a projection effect). This conclusion is reinforced by 



the emission-line spectra — the relative line fluxes of region a 
is almost identical to those of the other four regions (b-e), 
none of which are on the path of the radio jet (Table 3). 

As with 3C48, the [Oil] doublet is weU resolved in 
PKS 2251+1 1. But sufficient S/N is only reached in regions a 
and b. We measured [O ii] A3726/A3729 ratios of 0.97±0.08 
and 0.80±0.07 for a and b, respectively, which imply Ne = 
320±100 (380±120) and 130±70 (150±80) cm'^ if 7; = 10^ 
(1.5 xlO^) K. 

5. QUASAR PHOTOINIONIZATION AND GAS 
METALLICITY 

In Figure 9 we compare ratios of strong emission lines of 
the quasar EELRs with those of the emission-line galaxies in 
the Sloan Digital Sky Survey Data Release 4 (SDSS DR4). 
The emission-line fluxes, after correction for stellar absorp- 
tion and foreground Galactic extinction, are publicly avail- 
able for SDSS emission-line galaxies^. Following the same 
dereddening procedure as in Fu & Stockton (2008) and the 
classification schemes of Kewley et al. (2006b), we classi- 
fied the SDSS galaxies into Seyferts, LINERs, star-forming 
galaxies (H II), and AGN/star-forming composites. As can 
be seen in the BPT diagrams (Baldwin et al. 1981) shown 

^ http://www.mpa-garching.mpg.de/SDSS/DR4/; and see Tremonti et al. 
(2004) for a description of the data. 
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Fig. 10. — The same line-ratio diagrams as in Fig. 9, although ionization model grids have been plotted in replacement of SDSS Seyferts and LINERs. 
The dusty radiation-pressure dominated photoionization (Groves et al. 2004) model grids for gas-phase metallicities of 1/3 and 2/3 are shown in gray and 
black, respectively. Both grids assume a hydrogen density near the ionization front of 1000 cm~^, and they cover the same range of ionization parameters 
(0.0 > logf/ > -3.0, dashed lines with increasing thickness) and the same four power-law indices for the ionizing continuum {Ft, oc z/"; a = -1.2,-1.4,-1.7, 
and -2.0; solid lines with increasing thickness). The shock-only (blue) and "shock -i- precursor" models {red; Dopita & Sutherland 1996) show the same range 
of shock velocities (200 < V5 < 400 km s~^; dashed lines with increasing thickness) and four magnetic parameters {B/n^/'^ = 0.5, 1,2, and 4 fiG cm^/^; solid 
lines with increasing thickness). In panel b, the original photoionization model grids do not match the data. We applied a simple correction to the model grids 
by shifting them to the right by 0.51 in log, which is motivated by a recent update on the O/Ne abundance ratio [log (O/Ne) = 1.12 (Gloeckler & Geiss 2007) as 
opposed to the old value of 0.61 that was used in Groves et al. (2004)]. The girds after the correction are shown in thiner black and grey lines. 



in Fig. 9, the quasar EELRs show similar Hne-ratios to those 
of the Seyfert galaxies, but they are clearly distinguishable 
from star-forming and composite galaxies. This result, in 
combination with the inadequacy of shock and "shock + pre- 
cursor" models, as we demonstrated in Fig. 10, indicate that 
the EELRs are most definitely photoionized by an AGN-type 
ionizing spectrum. The only exception is the region 3C 48-a, 
which is below the AGN/star-forming dividing line in Fig. 9d 
and Fig. 11a. 3C48-a corresponds to a region between C5 
and C6 as defined in the long-slit observations of Canalizo & 
Stockton (2000a). C5 and C6 have roughly half of their stellar 
continuum coming from a young stellar population of only ~9 
Myr old. We believe that the low [Olll]/H/3 ratio of 3C48-a 
is due to a significant ionizing flux from the early-type stars 
in the region, i.e., 3C48-a is ionized not only by the central 
quasar > 26 kpc away but also by the local massive stars. 

Although we have compared the results of the dusty 
radiation-pressure dominated photoionization model (Groves 
et al. 2004) with our data in Fig. 10, we do not mean to 
imply that EELRs are actually dusty ionized surface layers 
of dense molecular clouds. We have demonstrated previ- 
ously (Fu & Stockton 2007b) that a two-phase photoioniza- 
tion model works equally well in reproducing the line-ratios. 



The dense ionization-bound clouds can easily be regenerated 
through turbulent shocks that should be prevalent in EELRs, 
given the observed supersonic velocity dispersions. Further- 
more, it is difficult to understand how dense molecular clouds 
could be present at large galactic radii, in the absence of any 
detectible underlying stellar continuum (Stockton et al. 2002). 

Some of the EELR clouds show noticeably weaker 
[S ll]/Ha ratios than those predicted by the photoionization 
models (Fig. lOd). [Six] (like [Ol]) depends strongly upon 
the density structure of the cloud, since it is mostly emitted by 
the partially ionized region behind the ionization front. There- 
fore, the observed inconsistency probably implies that the 
density structure of EELR clouds is different from a radiation- 
pressure dominated NLR cloud, offering further evidence 
against the dense molecular cloud model. Because of the un- 
certainties in modeling the density structure of the NLR cloud. 
Groves et al. (2006) have discounted [S ii] and [O i] ratios as 
reliable metallicity indicators. Instead, the authors concluded 
that "line ratios involving [Nil] ratios are the most robust 
metallicity indicators in galaxies where the primary source of 
ionization is from the active nucleus," mainly because of the 
secondary production of the element. On the other hand, by 
comparing two power-law photoionization models assuming 
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Fig. 11. — Quasar EELRs are mostly more metal-poor than the nebulae in typical AGN, as indicated by the metallicity-sensitive diagram of [O ill] X5001/¥lf3 
vs. [N II] A6583/Hq;. The background image shows the density distribution of the SDSS DR4 emission-line galaxies in log scale. Objects above the upper green 
curve are objects dominated by AGN (As in Fig. 9d, LINERs are concentrated in the lower denser branch and Seyferts in the upper branch; Kewley et al. 2006b), 
below the lower curve are star-forming galaxies (Kauffmann et al. 2003), and AGN/star-forming composite galaxies are in between. In panel (a) we show our 
measurements from the EELRs; while in panel (b) we show photoionization model grids (Groves et al. 2004) of four gas-phase metallicities from 0.2 to 1.3 
Zq. For each metallicity, model predications are given for a range of ionization parameters (-2.3 < \og(U) < 0) and four power law indices representing the 
quasar ionizing continuum (Fu (xiy^;a = -1.2,-1.4,-1.7, and -2.0). Most Seyferts can be fit quite well with the super-solar metallicity models. The top region 
bordered by the black dashed line and the upper green curve is where the low-metallicity Seyferts are located (Groves et al. 2006). 



drastically different density structures, we have shown that the 
[N II] Hne ratios are insensitive to the density profile but are 
almost entirely dependent on the nitrogen abundance (Fu & 
Stockton 2007b). 

In Figure 11 we estimate the metallicity of the gas in 
the EELRs with the metalHcity- sensitive [O III] A5007/H;^ vs. 
[Nil] A6583/Ha diagram. We have converted the modeled 
total metallicities to gas-phase metallicities using the solar 
abundances defined by Anders & Grevesse (1989), in order 
to be consistent both with previous studies on the metallic- 
ity of quasar BLRs and with metallicity results from stellar 
populations. One solar metallicity in Groves et al. (2004) cor- 
responds to 1/3 Zq in our figures, since approximately half 
of the metals are assumed to be depleted onto dust in their 
models, and their assumed solar metallicity is about 0.2 dex^^ 
lower than that of Anders & Grevesse (1989). In any case, in 
our comparisons with other AGN, it is the relative metallici- 
ties that are important, not the absolute values in solar terms. 

In agreement with our previous results (Fu & Stockton 
2008, their Fig. 12), all of the EELRs, except that of the only 
radio-quiet quasar, Mrkl014, are consistent with photoion- 
ization models with sub- solar metallicities (using the Anders 
& Grevesse 1989 scale), which are significantly lower than 
the metallicities deduced from typical AGN in the SDSS that 
span a similar redshift range. For Mrk 1014, metal-rich gas is 
detected in the western part of the EELR (regions a, h, and 
J), while the eastern cloud (c) shows a low metallicity similar 
to that of the other EELRs. Despite the precaution that [O I] 
and [S II] ratios may not be reliable metallicity indicators, we 
note that the EELRs also lie clearly to the left side of the 
AGN branch in the BPT [O l]and [S ll] diagrams (Fig. 9c and 
J), and they roughly overlap with the low-metallicity type-2 
Seyferts that were selected in SDSS by Groves et al. (2006). 

Recall that EELR quasars also distinguish themselves from 
other quasars by their BLR metallicities (Fu & Stockton 
2007a) — the former show much lower metallicities (Z < 

The C, N, O, and Fe abundances in the solar abundance set assumed by 
Groves et al. (2004) are all ~0.2 dex lower than those in Anders & Grevesse 
(1989). Hence, the mass fraction of metals (i.e. Z) is about 1.5 times lower 
— Zq= 0.013 — than that of Anders & Grevesse (1989). 



0.6 Zq) than the latter (> 1 Zq). The BLR metalHcities 
are estimated from broad permitted UV line ratios, such as 
N V A1240/C IV A 1549 and N v/He ii A 1640, based on the lo- 
cally optimally emitting cloud model (Baldwin et al. 1995). 
The BLR result is therefore based on completely different data 
and photoionization models than is the EELR metallicity, yet 
both indicate a low metallicity, if we adopt the most straight- 
forward interpretation for the low [N ll]/Ha ratios. Therefore, 
although neither the BLR model nor the NLR/EELR model is 
perfect, a lower gas metallicity of EELR quasars most consis- 
tently explains the totality of the data. 

6. SUMMARY AND CONCLUSIONS 

Here we bring together our results from IFU spectroscopy 
of the EELRs around 7 quasars and 1 FR II radio galaxy, as 
have been presented in this and previous papers (Fu & Stock- 
ton 2006, 2007b, 2008; Stockton et al. 2007). We summarize 
these results for several key properties and interpret them in 
terms of a tentative "toy model" for the formation of quasar 
EELRs. 

6. 1 . Physical Properties of Extended Emission-Line Regions 
6.1.1. Kinematics 

The velocity structure of the ionized gas is locally ordered 
but globally disordered in all of the 7 quasars, consistent with 
the results of previous work (Durret et al. 1994; Crawford & 
Vanderriest 2000; Stockton et al. 2002). We are not aware 
of any simple dynamical models that can fit the data. The 
only velocity map that resembles a rotating disk is that of the 
radio galaxy 3C79; however, even in this case, this is true 
only when the fainter clouds are ignored. 

The radial velocity maps have shown that most of the ex- 
tended nebular emission is confined within 300 km s~^ of the 
quasar's systemic velocity, which is defined by the nuclear 
NLR. Nevertheless, extremely high velocities are observed 
in fainter or presumably less massive clouds in 5 of the 7 
quasars— 3C 249.1 (+550 km s'^), 4C 37.43 (-600 km s'^), 
Mrk 1014 (-1100 km s'^), 3C48 (-500 km s'^), and Ton 616 
(+440 km s~^). We believe such high velocities are unlikely 
to be solely gravitational in origin. 
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The velocity dispersion maps look quite quiescent with typ- 
ical values between 30 and 100 km s~^ . Three of the 7 EELRs 
show clouds with velocity dispersions significantly greater 
than 100 km s"^: in 3C 249.1 there is a cloud to the SE of 
the nucleus with a ~ 310 km s~^ ; in Mrk 1014 there are two 
clouds to the E and W of the nucleus with cr ~ 280 km s~^ ; 
in 3C48 the high velocity cloud 0.''25 north of the nucleus 
shows a dispersion of ~390 km s~^ In all of the three cases, 
there seems to be a clear interaction between a radio jet and 
the cloud: the 3C 249.1 cloud is aligned with the radio jet 
to the SE, the Mrk 1014 clouds coincide with the two radio 
knots on either side of the core, and the high- velocity cloud 
in 3C 48 is near the base of the Q'.'l radio jet. However, no 
such kinematic signatures were seen in PKS 2251+11 where 
a positional correspondence between a radio hot spot and one 
of the brightest EELR cloud has long been noticed. Thus this 
is most likely a chance alignment; this conclusion is also sup- 
ported by the lack of evidence of shocks from line intensity 
ratios. 

6.1.2. Ionization Mechanism 

There is now little doubt that EELRs are photoionized by 
the central quasars. We reached this conclusion because pho- 
toionization by a power-law continuum source is the only 
model that fits the data. 

Shock ionization and the "shock + precursor" model cannot 
predict the correct line ratios; this result is especially clear in 
the [Oiii]/H^ vs. [Oii]/H^ diagram (Fu & Stockton 2007b, 
their Fig. 5a). The high [Olll]/H/3 ratios of the EELRs can 
only be reproduced by the "shock + precursor" model when 
the shock speed exceeds 450 km s~\ which is inconsistent 
with the low velocity dispersions measured in the clouds (typ- 
ically 30 < cr < 100 km s~^). Further, the model predicts an 
[O ll]/H/3 ratio ~ 2x higher than the observed values. 

Photoionization by star-forming regions has also been ruled 
out as a dominant ionization mechanism. For instance, the 
EELRs lie firmly above the maximum star-forming curve in 
the BPT diagrams of [O lll]/H/? vs. [S ll]/Ha and [O lll]/H^ 
vs. [Nii]/Ha (e.g.. Figs. 9 and 11). The cloud 3C48-a rep- 
resents the only case where in addition to the quasar young 
stars may also have contributed significantly to the ionization 
(§ 5). 

6.1.3. Metallicity 

EELRs consist of clouds with low metallicities (Z < 2/3 
Z0, on the Anders & Grevesse 1989 scale) as implied by 
their low [N ll]/Ha line ratios. These clouds are significantly 
more metal-poor than the gas in similar environments, e.g., 
the NLRs of type-2 Seyfert galaxies and radio galaxies at 
similar redshifts. Power-law photoionization models with low 
metallicities can correctly predict fluxes of optical lines other 
than [Nil] as well (see Fig. 10). Four of the seven quasars 
in our sample were observed by the HST EOS in the UV 
(3C 249.1, Ton 202, 4C 37.43, and PKS 2251+11), both the 
Nv A1240/CIV A 1549 and Nv/Heii A 1640 ratios of their 
BLRs indicate a metallicity less than 0.5 Z0(Fu & Stock- 
ton 2007a), remarkably consistent with the metallicity de- 
rived from the extended gas in a completely different density 
regime. 

The only radio-quiet quasar, Mrk 1014, is also the only one 
where extended metal-rich gas is detected. It seems that the 
gas in this quasar has not been well mixed, as a bright EELR 
cloud to the east shows a low metallicity that is similar to that 
of the other EELRs in the sample. We suspect the clouds to 



the west have been enriched by recent star-forming activities. 
Evidence for a recent starburst in the host galaxy of Mrk 1014 
was reported by Canalizo & Stockton (2000b). 

6.2. Star Formation in the Host Galaxies of EELR Quasars 

The host galaxies of a couple of our EELR quasars have 
been studied in great detail by Canalizo & Stockton (2000a,b). 
3C48 and Mrk 1014 were included in their sample because 
they fall into an intermediate region between quasars and 
ULIRGs in a FIR color-color diagram (Lipari 1994), suggest- 
ing significant star-forming components in their host galax- 
ies. Deep imaging of their host galaxies has shown clear 
evidence that these are cases of major mergers in their final 
stages. Deep spectroscopy revealed young stellar populations, 
confirming that FIR colors are indicative of star formation. 
In 3C48, Canalizo & Stockton found vigorous ongoing star 
formation and post- starburst populations with ages up to 114 
Myr. In Mrk 1014, only post- starburst populations with ages 
between 180 and 290 Myr were seen. Notice that 3C48 and 
Mrk 1014 are the only two objects in our sample that do not 
have classical FR II double-lobe morphologies. 

Similar studies do not exist for the rest of the quasars in 
our sample. Nevertheless, an additional two EELR quasars 
have been detected at 60 jim by IRAS and one by ISO — 
3C 249.1, 4C 37.43, and 3C 323.1. Although solid detec- 
tions are only available at 25 /im and 60 /im, their large 
a(60, 25) spectral indices place them firmly above the transi- 
tion region defined by Canalizo & Stockton (2001). The FIR 
colors^^ are a(60,25) = -0.35,-0.23,-0.05 and a(100,60) > 
-0.63,-2.10,-2.24 for 3C 249.1, 4C 37.43, and 3C 323.1, re- 
spectively. If ULIRG-unlike FIR colors are indicative of no 
recent star formation, then we would expect that the host 
galaxies of these EELR quasars will be devoid of young stel- 
lar populations less than a few Myr old. An indirect proof 
of this conjecture comes from our study of the host galaxy of 
3C79 (Fu & Stockton 2008). 

Given that 3C 79 is a type-2 EELR quasar with FR II mor- 
phology, the unique geometry of the radio galaxy offers a 
much clearer view of the host galaxy of an EELR quasar. We 
found that the host galaxy consists of an intermediate- age (1.3 
Gyr) stellar population (4% by mass) superimposed on a 10 
Gyr old population, i.e., there has been no recent star forma- 
tion. 

In summary, the host galaxies of EELR quasars can have 
a large range of star formation histories. For Mrk 1014 and 
3C 48, stars have been vigorously forming in the past few hun- 
dred Myr; but for 3C 79, and, apparently, for at least many of 
the FR II quasars in our sample, significant star formation has 
occurred only in the remote past (> 1 Gyr ago). 

6.3. A Toy Model of Quasar Extended Emission-Line 
Regions 

A scenario for producing quasar EELRs must be able to ac- 
count for all of the following compelling results: (1) the corre- 
lation between extended nebular emission and radio spectral 
index (§ 2.1), (2) the correlation between extended nebular 
emission and quasar BLR metallicity (Fu & Stockton 2007a), 
(3) the lack of detailed morphological correspondence be- 
tween extended nebula and the host galaxy or the radio struc- 
ture, (4) the chaotic velocity distribution but relatively calm 

When calculating these spectral indices, photometric data from both 
IRAS and ISO have been considered. At 25 iim and 60 /im the priority list is 
IRAS soHd detection > ISO solid detection > IRAS upper limits, and at 100 
lim only IRAS upper limits were used. 
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velocity dispersions of the extended nebulae, (5) the frequent 
detection of extremely high velocity clouds (V > 400 km s~^), 
and (6) the large masses of the EELRs (M - 10^-10^^ M© ; Fu 
& Stockton 2006, 2007b, 2008; Stockton et al. 2007). 

As mentioned in § 2.2, there were three remaining at least 
somewhat plausible origins for the extended gas, and we now 
examine each one of them with the current set of observational 
results: 

1. Cold accretion of intergalactic gas through Mpc-scale 
filaments. Although this model can easily explain the 
low metallicity gas comprising the EELRs, it is difficult 
for it to explain the quasar activity and the similarly 
low metallicity of the BLRs. Cold accretion by itself 
is unlikely to activate the central black hole as there is 
no likely mechanism to provide the drastic reduction of 
the angular momentum of the accreted gas that would 
be required to force a significant amount of it to within 
the accretion radius of the black hole. 

2. Tidal debris from a galactic merger. It is difficult for 
this model to explain the clouds with velocities higher 
than reasonable for gravitational kinematics (> 400 
km s~^), which have been detected in 5 of the 7 quasars, 
often at a considerable distance from the quasar and 
in locations not in the path of a radio jet (3C 249.1, 
Ton 616, 4C 37.43). In addition, it is difficult to see 
how this tidal picture would give a situation where 
no luminous EELRs have been found associated with 
flat-spectrum core-dominated quasars, and those asso- 
ciated with radio-quiet QSOs are relatively rare and 
of generally lower luminosity than those found around 
steep- spectrum radio sources. Nevertheless, most of the 
quasars with luminous EELRs do show independent ev- 
idence for recent strong interactions or mergers (Stock- 
ton & MacKenty 1987), and this evidence has to be 
taken into account in any reasonable physical model. 
But it is significant that we never see any indication 
of an underlying stellar component with a distribution 
similar to that of a luminous EELR. 

3. Remnants of galactic superwinds driven by either the 
quasar itself or a coeval starburst. The starburst su- 
perwind picture has been found problematic both in 
3C 249.1 on the ground of the derived momentum in- 
jection rate (Fu & Stockton 2006) and in 3C 79 based 
on the lack of young stellar populations (Fu & Stockton 
2008). In addition, starburst-driven superwinds reflect 
mass ejection from type-II supernovae, and thus they 
are expected to be metal enriched (it takes only ~0. 1 
Gyr to reach a gas metallicity of 1 Zq from the primor- 
dial abundances; Friaca & Terlevich 1998). This con- 
tradicts the dominance of metal-poor gas in the environ- 
ments of EELR quasars. It would be especially difficult 
to explain the metal-poor gas in the quasar BLRs, which 
are within 0.1 pc of the central black holes. Finally, a 
starburst-driven superwind cannot explain the fact that 
luminous EELRs are found virtually exclusively around 
steep- spectrum radio-loud quasars. On the other hand, 
EELRs as remnants of superwinds driven by quasars 
themselves do not suffer from these difficulties. 

We clearly need to combine aspects of the two last points 
into our over-all picture. Both strong interactions or mergers, 
on the one hand, and the presence of a steep- spectrum radio 



source, on the other, each in general appear to be necessary, 
but not sufficient, ingredients for the production of a luminous 
EELR. There is considerable evidence that the merging com- 
panion supplies the gas (Fu & Stockton 2007a), but that the 
distribution of the ionized gas is controlled by a massive out- 
flow (Stockton et al. 2002, 2007; Fu & Stockton 2006, 2007b, 
2008). We suggest the following overall picture for the forma- 
tion and the fate of the EELRs around quasars with large FR 
II radio jets: A large, late-type galaxy with a mass of a few 
X 10^^ M0 of low-metallicity gas has recently merged with 
the gas-poor quasar host galaxy, which has a ~ 10^ Mq black 
hole at its center. Low metallicity gas on the outskirts of the 
merging companion is driven to its center and is completely 
mixed with any higher metallicity gas there well before the 
final coalescence (e.g., Kewley et al. 2006a). The merger trig- 
gers the current episode of quasar activity, including the pro- 
duction of FR II radio jets. The initiation of the jets also pro- 
duces a wide-solid-angle blast wave that sweeps most of the 
gas from the encounter out of the galaxy. This gas is imme- 
diately photoionized by UV radiation from the quasar when 
the blast wave clears the dust cocoon enshrouding the central 
engine. Turbulent shocks produce high-density (~ 400 cm~^) 
filaments or sheets in the otherwise low-density (~ 1 cm~^) 
ionized medium which is in hydrodynamical equilibrium with 
a hot (T ^ 10^ K) diffuse medium, maintaining a low ioniza- 
tion state. The quasar will eventually shut itself off due to 
starvation, and some high velocity clouds will escape into the 
inter-galactic medium, but most of the warm gas will stay and 
recombine to neutral hydrogen (HI) ~ 10"^ yr after the quasar 
shuts down. The chaotic kinematics will be washed away as 
the gas settles into a rotating disk in a crossing time of ^0.5 
Gyr. But before it settles down, it might look much like the 
HI gas in the nearby early-type galaxy NGC 1023 (Morganti 
et al. 2006). The dynamically settled gas could probably ac- 
count for a significant fraction of the dispersed HI in giant 
early-type galaxies. 

Note that, in 3C48, we may be witnessing an early phase 
in one of these episodes (§ 4.1; Stockton et al. 2007). The 
high-velocity gas currently concentrated near the base of the 
radio jet constitutes a wide-solid-angle outflow with velocities 
ranging up to at least 800 km s~^ and a total mass of ionized 
gas of at least ~ 10^ M©, and quite likely much more. 

This picture naturally explains the EELR — radio morphol- 
ogy correlation, since the superwind has been assumed to be 
associated with the formation of the jets. But the reason for 
the EELR — low gas metallicity correlation is much less clear. 
We can speculate that it may have something to do with the 
lower radiative coupling such gas would have to the quasar ra- 
diation field, allowing more efficient accretion. Specifically, 
a higher metallicity will lower the accretion rate of material 
towards the center, because both the higher opacity of the gas 
and larger amount of dust will couple the gas more efficiently 
to the radiation field of the quasar. Such a lowered accretion 
rate may delay the spin up of the black hole, consequently de- 
laying the formation of the radio jet, assuming that the jet pro- 
duction is determined by the black hole spin (Blandford et al. 
1990). Hence, most of the gas may have time to form stars be- 
fore the jet is launched, leaving an insufficient amount of gas 
to form an EELR. It is interesting to note that related consider- 
ations have been invoked recently to explain the link between 
low metallicity and long duration 7-ray bursts (Fruchter et al. 
2006). 

Whether this particular explanation for the low metallici- 
ties in the BLRs of EELR quasars is correct or not, the cor- 
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relation itself clearly suggests some physical mechanism by 
which the injection of a large amount of low-metallicity gas 
into a system with a super-massive black hole can result in a 
large-scale, quasar-powered superwind. This sort of situation 
is likely to have been relatively common in the early universe, 
when mergers were more frequent and gas-rich galaxies more 
numerous. This form of feedback might therefore have been 
of some importance in restricting the growth of the most mas- 
sive galaxies at early epochs. 
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TABLES 

Line Ratios of Quasar Emission-Line Clouds Relative to }if3 



ID 


Ay 


[Nev] A3426 


[O II] A3727 


[Neiii] A3869 


[Oiii] A4363 


HeiiA4686 


[Oiii] A5007 


[Oi] A6300 


Ha 


[Nil] A6583 


[S II] A6716 


[S II] A6731 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


3C 48 (z = 0.369) 


a 


(2.80) 




4.63±0.19 


0.67 ±0.14 




1.00±0.04 


3.63 ±0.04 


0.12 ±0.02 


3.10±0.08 


0.57 ±0.06 


0.35 ±0.06 




b 


(2.80) 




6.25 ±0.27 


1.94 ±0.21 




0.42 ±0.03 1.00 ±0.03 


6.67 ±0.02 


0.14 ±0.02 


3.10±0.10 


0.50 ±0.08 


0.31 ±0.06 




c 


2.80 


2.80±0.18 


4.63±0.11 


1.56 ±0.05 


0.21 ±0.02 


0.43±0.01 1.00±0.01 


7.58 ±0.01 


0.16±0.01 


3.10±0.03 


0.39 ±0.02 


0.23 ±0.03 




d 


(2.80) 


1.54±0.17 


4.63±0.13 


1.40 ±0.05 


0.21 ±0.02 


0.42 ±0.01 1.00 ±0.02 


8.64 ±0.01 


0.19±0.01 


3.10±0.03 


0.50 ±0.02 


0.25 ±0.02 




Mrk 1014 (z = 0.163) 


a 


0.61 




3.10±0.21 


0.83±0.11 




0.46 ±0.03 1.00 ±0.04 


7.61 ±0.03 


0.25 ±0.02 


3.10±0.02 


2.38 ±0.02 


0.65 ±0.01 


0.50 ±0.01 


h 


0.34 




2.93 ±0.25 


1.04 ±0.12 




0.38 ±0.03 1.00 ±0.04 


8.06 ±0.04 


0.32 ±0.02 


3.10±0.02 


2.91 ±0.02 


0.80 ±0.01 


0.62 ±0.01 


c 


(0.00) 


1.48 ±0.34 


1.98±0.17 


0.55 ±0.06 




0.47 ±0.02 1.00 ±0.02 


9.07 ±0.03 


0.13±0.01 


2.86 ±0.01 


1.04 ±0.01 


0.47 ±0.01 


0.36 ±0.01 


d 


1.15 




3.94 ±0.35 


1.48 ±0.26 




0.42 ±0.06 1.00 ±0.07 


4.64 ±0.06 


0.26 ±0.02 


3.10±0.02 


2.65 ±0.02 


0.70 ±0.03 


0.46 ±0.02 


Ton 616 a = 0.268) 


a 


0.81 










0.56 ±0.03 1.00 ±0.03 


11.05 ±0.03 




3.10±0.03 


0.70 ±0.02 


0.15 ±0.02 


0.14 ±0.02 


b 


0.59 










0.31 ±0.03 1.00±0.03 


10.14 ±0.02 


0.09 ±0.02 


3.10±0.02 


0.61 ±0.02 


0.44 ±0.02 




c 


0.28 










0.42 ±0.09 1.00 ±0.10 


11.99 ±0.07 


0.26 ±0.03 


3.10±0.05 


0.58 ±0.04 


0.39 ±0.08 




d 


0.14 










1.00±0.07 


12.68 ±0.05 


0.16 ±0.02 


3.10±0.03 


0.52 ±0.02 


0.31 ±0.03 




e 


(0.00) 










0.26 ±0.05 1.00 ±0.06 


10.81 ±0.05 


0.19 ±0.02 


3.00 ±0.03 


0.79 ±0.02 


0.48 ±0.03 


0.36 ±0.03 


Ton 202 (z = 0.364) 


a 


(2.30) 


1.34 ±0.20 


3.02±0.16 


0.94 ±0.12 




0.37 ±0.03 1.00 ±0.03 


7.35 ±0.02 




3.10±0.07 


1.11±0.06 


0.28 ±0.04 


0.27 ±0.04 


b 


(2.30) 


1.59 ±0.29 


3.01 ±0.30 


0.75 ±0.12 




0.42 ±0.05 1.00 ±0.05 


7.22 ±0.04 




3.10±0.15 


0.53±0.12 


0.38 ±0.09 




c 


(2.30) 


2.28 ±0.56 


4.66 ±0.56 


1.38 ±0.25 




0.43 ±0.07 1.00 ±0.07 


11.41 ±0.07 




3.10±0.22 




0.93 ±0.20 




PKS 2251+11 (z = 0.325) 


a 


(0.46) 




2.94 ±0.12 


0.97 ±0.07 




0.37 ±0.03 1.00 ±0.04 


10.49 ±0.04 




3.10±0.05 


1.10±0.04 


0.51 ±0.05 




b 


(0.46) 


1.01 ±0.21 


3.52±0.14 


1.08±0.12 




1.00±0.05 


9.12 ±0.04 




3.10±0.06 


1.49 ±0.05 


0.60 ±0.06 


0.49 ±0.05 


c 


(0.46) 




4.69 ±0.17 


0.88±0.11 




0.21 ±0.05 1.00 ±0.06 


8.79 ±0.03 




3.10±0.06 


1.87 ±0.05 


0.70 ±0.05 


0.49 ±0.04 


d 


(0.46) 


0.89 ±0.21 


2.84 ±0.23 


0.44 ±0.09 




0.26 ±0.05 1.00 ±0.06 


7.16 ±0.04 




3.10±0.07 


1.41 ±0.06 


0.50 ±0.08 




e 


(0.46) 


1.32 ±0.30 


2.78 ±0.41 


0.86±0.13 




1.00±0.07 


9.99 ±0.06 




3.10±0.06 


1.47 ±0.06 


0.59±0.12 


0.50±0.11 



Note. — Col. (1): Cloud name. Col. (2): Intrinsic reddening. The values in parentheses are not directly derived from the Balmer decrements of the corresponding clouds. The Ay values are measured from the Hq!/H/3 ratio 
for Mrk 1014 and Ton 616, where the two lines were observed simultaneously. For 3C48, Ton 202 and PKS 2251+11, they have been fixed respectively to those derived from the H7/H/3 ratios of clouds 3C48-C, Ton202-a + /?, and 
PKS 2251+1 1-a + ^ + c, because the S/N ratio of an individual cloud at H7 is normally insufficient for a meaningful measurement of the Balmer decrement. Ay with derived values slightly less than zero have been fixed to zero. Col. 
(3-14): Dereddened intensities of emission lines relative to H/3. The line-of-sight Galactic extinction has also been taken into account. Undetected lines (S/N < 4) have been omitted. 



